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Abstract 

In a recent paper we have presented a method to evaluate certain Hankel determinants as 
almost products; i.e. as a sum of a small number of products. The technique to find the explicit 
form of the almost product relies on differential-convolution equations and trace calculations. In 
the trace calculations a number of intermediate nonlinear terms involving determinants occur, 
but only to cancel out in the end. 

In this paper, we introduce a class of multilinear operators 7 acting on tuples of matrices 
as an alternative to the trace method. These operators do not produce extraneous nonlinear 
terms, and can be combined easily with differentiation. 

The paper is self contained. An example of an almost product evaluation using 7-operators 
is worked out in detail and tables of the 7-operator values on various forms of matrices are 
provided. We also present an explicit evaluation of a new class of Hankel determinants and 
conjectures. 

Mathematics Subject Classifications: 05A10, 05A15, 05A19, 05E35, 11C20, 11B65 



1 Introduction 

The expansion of a determinant 

det [aij]o<ij<n 

from first principles involves calculating the signed sum of (n + 1)! individual products. This type 
of an evaluation is not of much interest, and one usually uses the multilinearity of the determinant 
to obtain more succinct expressions for a given family of determinants. Those determinants which 
may be evaluated as a single product of simple factors (such as the Vandermonde and Cauchy 
determinants) have a special appeal. For product form evaluations, LU decomposition, continued 
fractions and Dodgson condensation are some of the available methods that have been utilized with 
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considerable success. There exists an extensive literature on this topic, going back to the treatise 
of Muir [lOl [11] . A more recent compilation of the state of affairs of the theory of determinants 
appears in Krattenthaler [7] , in which a wide range of techniques used to study the evaluation of 
families of determinants are described, accompanied by an extensive bibliography on the subject. 
Of particular interest are Hankel determinants, for which 



H+j 



Certain classes of Hankel determinants with combinatorially interesting entries ai^j have product 
representations with startling evaluations, and we mention 



det 



and 



det 



Hi + j) + 2 
i+j 



3(»+j) 
i+j 



n 



(6i + 4)!(2i + 1)! 
o<.,<„ t\ 2(47T2)!(47+3)! 

3(3i + l)(6i)!(2i)! 



n 



(4i)!(4i + 1)! 



0<i.j<n i—\ 

(see [1] and [3]). The in-between case of the binomial coefficients 



3fc + 1 
k 



(1) 



is not amenable to standard methods since it does not have a product evaluation. In a recent paper 
[2] we proved that for the entries ((T]), the evaluation is an almost product] in this case a sum of n + 1 
products of simple factors: 



det 



3(i + j) + l 
i+j 



^ (6i + 4)!(2i + 1)! ^ 7i!(4n + 3)!!(3n + i + 2)! 
o<..j-<n ~ tX 2(4* + 2)!(4z + 3)! ^ (3n + 2)H!(,i - z)!(4n + 2z + 3)!! 



The technique presented in [5] to find the explicit form of the almost product for this particular 
Hankel determinant relies on the following steps: 



(I) Using k = i + j, replace with polynomials 

k 



ak 



3fc + 1 — TO 

k — m 



(2) 



so that afc(x) is a monic polynomial of degree k with Ofe = 0^(0). Consequently the associated 
Hankel determinant Hn{x) is a polynomial, and = Hn{0). 

(H) Establish a second order ODE satisfied by Hn{x). 

(HI) Solve the DE in (II), and evaluate the solution at a; = 0. 

The {P,a)-case of this problem is the evaluation of the Hankel determinants where the entries 



are 



-,(/3,") 



m=0 



f 13k + a — m 
\ k — m 



(3) 
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The bulk of the work is contained in Step (II) , and this part of the argument itself relies on three 
essential identities. These identities are linked in the derivation of the differential equation via the 
application of a trace operator. 

In this paper, we introduce a class of multilinear 7-operators acting on tuples of matrices which 
take the place of this trace operator. 

If it had just been a matter of calculating the differential equation in the (3, l)-case as we did 
in [2], then which technique we used might not have mattered much. However, we wanted to try to 
extend the differential equation method to a larger class of (/?, a)-cases, and we found that already 
in the (2, 2)-case, the 7-operators simplified the calculations significantly. To be specific, in the 
trace approach some nonlinear terms occur in the calculations, which get canceled in the end. For 
example the following ratio of determinants (using the notation in [2]) 



appears during the course of the trace calculations (e.g. p;, p. 15), and is later cancelled. 

As one goes to other cases, these nonlinear terms proliferate. In the (2, 2)-case, there are over 
half a dozen of these terms that arise, which all cancel. 

These nonlinear terms turn out to be an avoidable burden in a method that already involves a 
lot of calculation. It is easier to combine differentiation with the 7-operators than with the trace cal- 
culations of [2] and in addition the 7-operator calculations do not produce the extraneous nonlinear 
terms mentioned above. An added benefit is that they need not be calculated from scratch for other 
Hankel determinant evaluations. In Appendix III, we provide extensive table of values of 7-operators. 



4(4n-f 3)2-Ii 



(4) 



Let 




k 



(5) 



and define the {n + 1) x {n + 1) Hankel determinants by 



H,i{x) = det[ai+j{x)]o<i,j<n 



(6) 



A few of these polynomials and the Hankel determinants are as follows: 



ao{x) 1 

ai{x) — 4 + a; 

o-iix) = 15 + 5a: + 

a-i{x) = 56 + 21x -I- Bx^ -I- 



and 



H^{x) = 1 

Hx{x) = -1-32: 

H2{x) = -\-x + hx^ 

H-i{x) = 1 + 6x -I- 3x2 - . 
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We give the elements of the appUcation of 7-operators by working through the proof of the following 
theorem. 



Theorem 1 Suppose and the Hn{x) are as defined in and (S^. Then Hn{x) has the following 
almost product evaluations: 



k=0 



and 

n 



n + k\ , ^ fn + k + 1 



{X-2Y 



k=0 



{n + k + l)[ +(2n + 4fc + l)L^J+8(fc + l) ^ ^ 



(7) 



{x+2 f . 
(8) 



Alternate expressions for ([7]) and ([8]) are given in ((43l) and (|46)) . The expansion of Hn{x) around 
X = can be found in ([51]) . The generating function of the iJ„ (x) itself is given in (|44p . 
It is known that [1 [H] 



det 



2{t+j) + 2 
i+j 



(-1)- 



(9) 



0<^,J<n 



Our purpose is not the derivation of this relatively simple numerical evaluation itself, but to give 
an exposition of the salient points of the 7-operators, which allow us to evaluate the general case of 
the Hankel determinants of the polynomials ([5]) as an almost product. 

Additionally, we obtain numerical evaluations of Hn{x) at special values of x. A number of these 
are presented in Section [5] and at the end of Section [51 

In Corollary [3] we evaluate the Hankel determinant 



det 



2{i + + 
i + j 



0<«,j<n 

The explicit almost product evaluation of Theorem [T] is derived from the second order differential 
equation satisfied by these Hankel determinants. This differential equation is given in Theorem [2] 
in Section [71 With the definition of the polynomials in the evaluation in this paper is the 
(/?, a) = (2, 2)-case. 

The outline of the rest of this paper is as follows: In Section [21 we define determinants H\ for 
partitions A obtained from a given Hankel matrix. This is followed by the introduction of the family 
of multilinear operators 7 along with their basic properties and a combinatorial interpretation for 
their evaluation in Section[21 Section[31presents example calculations with the 7's, and a compilation 
of evaluations that are used in the paper. This is followed by three identities that are typically needed 
for our methods, and the derivation of the equations satisfied by the various H\ that arise in the 
calculations. We obtain a system of first order differential equations which results in a second order 
differential equation for the Hankel determinant we wish to evaluate in Section [71 Evaluation at 
special points are discussed in Section \8\ and the general solution of the differential equation is 
derived in Section [9l An additional Hankel determinant evaluation is given at the end of this section 
in Corollary [3l In Section [TOl we consider the properties of the zeros of the Hankel determinants and 
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show that they form a Sturm sequence. Conjectures on the evaluation of similar Hankel determinants 
are presented in Section [TlJ This is followed by Appendix I - III where we give the proofs of the 
results stated and used in the calculations as well as tables of 7-operator evaluations. We remark 
that Sections 2-4 and Appendices I and III apply to general Hankel matrices, whereas Sections 5-10 
and Appendix II apply to the evaluation of the case a = (5 — 2. 



2 Preliminaries 



We consider general Hankel matrices A — [ai+j\o<ij<n in the symbols afc. In [5] and in Section [T] 
of the present paper we used the notation iJ„ for det(A). However, it is useful to have alternate 
notation for various determinants that arise, in which sometimes the parameter n is suppressed. 
Unless otherwise indicated, we assume that n has been chosen and is fixed. 

A partition A of an integer m > is a weakly decreasing sequence of integers A ~ (Ai > A2 > 
■ • • > Ap > 0) with TO = Ai + A2 + • • • + Ap. Each Ai is called a part of A. For example A — (3, 2, 2) 
is a partition of m = 7 into p = 3 parts. 

We use the notation A = to"'" • • • 2"^ 1"^ for a partition A of m, indicating that A has ai parts of 
size i. Thus for example, A — 3^21^ denotes the partition 3 + 3 + 2 + 1 + 1 + 1 of 11. Given n > 0, 
each partition (Ai > A2 > ■ • • > Ap > 0) with p < n + 1 defines a determinant of a matrix obtained 
from the (n + 1) x (n + 1) Hankel matrix An = [ai+j\o<i,j<n in the symbols a^, by shifting the 
column indices of the entries up according to A as follows: Let /i; = Ai for i ~ 1, . . . ,p and /ii = 
for i=p+l,...,n+l. Then 

H\ = det[ai+j+^^^^_ Jo<i j'<n ■ 

We use the special notation to denote the sequence m — Q lor i — 1, . . . , n + 1. For example when 
n — i. 



Ho = det 



ao ai a2 03 

ai a2 03 04 

a2 03 04 as 

az 04 05 ae 



H2 = det 



flo fli a2 as 

ai a2 as og 

a2 as a4 ay 

03 04 as ag 



det 



ao a2 
ai as 
a2 a4 



as ae 
a4 ay 
as as 



a3 as ag aio 



We note that these determinants are obtained in a way similar to the expansion of Schur functions 
in terms of the homogeneous symmetric functions by the Jacobi-Trudi identity . 

When the a^ = afc(x) are functions of x, then Hx = Hx{x) is a function of x. When we need to 
indicate the dependence of the determinant on n as well as x, we write 

Hx{n,x) 

for the (n+ 1) X (n+ 1) shifted Hankel determinant. As an example, with this notation ([7|) is written 

2n2 + 471 + 2fc2 + 1 



as 



Ho{x)=Ho{n,x) - (-1)" 



k=0 



2k + 1 



71 + fc 

2k 



yx-2r 



(10) 



The (n + 1) X (tt, + 1) Hankel determinant will be denoted by a number of different notations in 
this paper. Among these are iJ„ — iJ„(x), Hq — Hq[x), and Ho{n,x). In the latter two cases it 
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should be clear from the context that the subscript refers to the partition involved and not to the 
dimension of the Hankel matrix. 

Our aim is to obtain a first order linear system of equations 

Q-^Ho = QoHo + QiHi (11) 

U-^Hi - UoHo + UiHi 
ax 

where the coefficients are polynomial functions of x and n. From this system the second order 
differential equation for Hq given in Theorem [5] can be found immediately. 
In the process of differentiating Hq and Hi the following five determinants 

are encoimtered. We will express each of these in terms of the two determinants Hq,Hi. 

The 7-operator that we next define allows us to do this from the three identities satisfied by the 
Uk, while avoiding having to deal with nonlinear expressions involving determinants. This operator 
has the additional advantage of simplifying differentiation of determinants, improving on the trace 
calculations used in [5]. 

3 The 7-operator 

We define a multilinear operator 7 on m-tuples of matrices as follows: 
Definition 1 Given {n + 1) x (n + I) matrices A and Xi, X2, ■ ■ ■ , Xm with m > 1, define 

7a( ) = det(A) 

and 

ja{Xi, X„,) = dt.dt, ■ ■ ■ dt,„ det{A + tiXi + t2^2 + • • • + t„,X„^)\t,=-=t,„=o 
where ii, ^2, ■ ■ • , are variables that do not appear in A or Xi, X2, . • . , X^- 

Next we give a computationally feasible combinatorial interpretation of 7a(^i, . . . , Xm) for small 
m, based on elementary properties of determinants. 

Definition 2 Suppose A and Xi, . . . , X„i are (n + 1) x (ri + 1) matrices, m < n + 1. Given a subset 
of column indices S = {ji, j2, ■ ■ ■ , jm} ^ {0, 1, . . . , n} and a permutation a of {1, 2, . . . , to}, As^^ is 
defined as the matrix which is obtained from A by replacing A's jk-th column by the jk-th column of 
the matrix X^^, for fc = 1, 2, . . . , to. 

With this notation we have 

Proposition 1 For to < n + 1, 

7a(Xi,...,X„) = ^det(As,.) (12) 

where the summation is over all subsets S of {0, 1, . . . ,n} with \S\ = to and all permutation a of 
{1,2,...,to}. 
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Note: The expansion (|12p is also valid as a sum over row indices where the replacements made are 
rows from Xi, . . . , Xm instead of columns. 

Another motivation for using the 7-operators is that they differentiate nicely; the derivative of a 
7 is a sum of 7s. 

Proposition 2 For m < n, 

d d d 

^7a(-'«^1, ■ • ■ ,Xm) = jAi-^A,Xl, . . . , X„i) + ^ jAiXi, . . . —Xj,Xj + i, . . . ,X„i) . 

The proofs of Proposition [T] and Proposition [5] can be found in Appendix I. 

Using Proposition [H we can evaluate 7^ on matrices that are associated with A in terms of 
determinants Hx for various partitions A. Next, we give a few examples of these calculations and a 
compilation of the expansions needed. 

4 Explicit 7^ evaluations 

Let 

We start with a few sample calculations. 

Example: In the calculation of 'yA{[ai+j]), the sum in (I12|) is over all subsets S C {0, 1, . . .71} 
with a single element and a is the identity permutation. We are replacing a column of A with the 
same column, so the resulting determinant is Hq = dot (A) for each one of n + 1 possible column 
selections. Thus 

7A(h+j]) = (" + l)Fo • 

Example: In the calculation of jA{[0'i+j+2]) the sum in is again over all subsets S C 
{0, 1, . . .n} with one element. If 5 = {j} and j < n — 2, then the j-th and the (j + 2)-nd columns are 
identical in As^a and the determinant vanishes. For j — n, the determinant is H2 and for j = n — 1 
it is —if 12. Therefore 

lA{[ai+j+2]) = -ff 2 - Hi2 . 
Example: We split the calculation of jA{[{i + j)a,i+j+2]) into two pieces: 

7A([(i +i)aj+j+2]) = JA{[ia^+J+2]) + JA{[ja^+J+2])■ 

In the calculation of 7_A([jcii+j+2]), the determinant in (|12p survives only for S — {n} and S = {n— 1}, 
exactly as in the case of the evaluation of 7A([ai+j+2]) above. However, now the determinant gets 
multiplied by the factor n of the new n-th column in the former case, and by the factor n — 1 of the 
(n — l)-st column in the latter. Therefore 

7A([ja,;+j+2]) = nH2 - {n - l)Hi2 
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jA{[icii+j+2]) evaluates to the same expression, since now we are dealing with rows instead of 
columns, but otherwise the argument is the same. Therefore 

jA{[{i + j)a;+j+2]) = 2nH2 - 2{n - l)Hi2 . 

Definition 3 For a polynomial sequence a„ = a„(a;) (n > 0), the convolution polynomials Cn = 
Cn{x) are defined by 

n 

— ^ ^ ^k^n — k 
k=0 

with C_i = 0. 



Example: To compute ^A{[ci+j+i]) for n = 2, we use the expansion of the matrix [Q+j+i] in terms 
of shifted versions of A as given below. The expansion for arbitrary n can be found in Appendix I. 





ai 


02 


03 






ao 


ai 


a2 









ao 


ai 










ao 




0.2 


as 


04 


+ 


ai 


ai 


a2 


as 


+ 







ai 


a2 


+ as 








ai 




. ^3 


04 


as 






a2 


as 


04 









a2 


as _ 




_ 





02 








































+ 


flo 


a2 


03 


04 


+ 


ai 



























. "3 


04 


as 






"2 


03 


a4 















A routine calculation gives 



lA{[ci+j+i]o<i,3<n) = aoHi + naiHo + a^Hi + naiHo 
= 2ao-ffi + 2naiHo . 

We provide another example of a 7 calculation. 
Example: In the calculation of 7A([ai+j+i], [ai+j+2]) the sum in is over all subsets 5* C 
{0, 1, . . . n} with two elements. If S' = {ji < ^2} with j 2 < ^ 2, then for a = (1)(2), the columns 
j2 and j2 + 2, and for <t — (12), the columns j2 and j2 + 1 of As\a are identical. Therefore in these 
cases the determinant vanishes. The remaining possibilities for S, a pairs can be enumerated as 

1. 5 = {n-l,n} and (7= (1)(2), 

2. 5 = {n- 2,n- 1} and a = (1)(2) , 

3. 5 = {n - 2, n} and a = (12) . 
The resulting determinants are 

H2I, ~ Hl3, — i?l3, 

respectively. Therefore 

7A([aj+j+i], [ai+j+2]) = -^21 - 2iJi3 . 

In Tables [211S] of Appendix III, we give a list of various 7 evaluations. The ones that are needed 
for the computations in this paper are in Tables [2] and [3l 
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5 The three identities 



Now we consider the (2, 2)-case. The three identities used in the argument are given in the following 
three lemmas. These identities are typical of our methods. The first identity is a differential- 
convolution equation. The second identity involves convolutions and ak but no derivatives. The 
third identity is a linear dependence among certain column vectors involving the a^. 



Lemma 1 (First Identity (FI)) 



d 



-{ix^ + 18a;2 - 20.x + 24 + 4n(.x2 + 4))a„ (14) 
+8(a; - ifcn - i2{x - ifcn-i 



(x - 2)a;(a; + 2)(3x + 2)— a„ = 2n(a; - l)a„+2 + - 6)(a; - 2) + 3a;^ - 2a; + 4)a„+i 
dx 



Lemma 2 (Second Identity (SI)) 



{nx + 3x + 2)a„+2 — {nx{x + 6) + ix^ + IGx + 8)a„+i 

+2x{x + 2){2n + 5)a„ + {x - l){x - 2)c„ - 4(a; - 1)(.t - 2)c„_i = (15) 
Lemma 3 (Third Identity (TI) ) 

n+2 

Wn,] {x)ai+j (x) = (16) 

j=0 



for i — 0,1, . . . ,n where 



r 2(2n + 5) fn + i + 2\ (2n + 3)(2n + 5) fn + j + 2 
I 2.7 + 1 2j j+ 27Tl V 2j 

(2n + 3)(2n + 5) (n + j + 2 



(17) 



2j + 3 V 2j + l 

The proofs can be found in Appendix II. We remark that the weights in Lemma |3] are typical of 
our method. Once the coefficients of the weight polynomials Wn,j{x) are guessed, then automatic 
binomial identity provers can be used to prove ([T| 



To prove pT|) . we will find the expansions of both -^Hq and -^Hi in terms of Hq and Hi. Since 
at first other determinants II\ also appear in these derivatives, they will need to be eliminated. We 
do this by constructing a sufficient number of equations involving them, and then expressing each 
one in terms of Hq and Hi . 

6 The five equations 

6.1 Equation from 7^([S'/(i + j)]) 
Apply 

1a{*) 
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to the (n+1) X (n+1) matrix whose (i,j)-th entry is obtained from the second identity (|15p evaluated 
at i + j and expand using hnearity. If we denote the matrix so obtained from the second identity 
by [SI{i + j)], then the computation is the expansion of 7yi([5'/(i + j)]) = 0. We obtain 

-x{x + 6)jA{[{i + j)at+j+i]) - (3x^ + 16x + 8)jAi[at+j+i]) 
+4x{x + 2)jA{[{i + j)a^+J]) + mx{x + 2)JA{[a^+J]) 
+ {x -l){x~ 2)7a([q+j]) - 4(x - l)(a; - 2)7a([c,+j_i]) 
Making use of the entries in the 7_a(*) computations from Table [21 we get 
= x{2nH2-2{n-l)Hi2) + {3x + 2){H2- H12) 
-x{x + Q)2nHi - (3x2 ^ -^g^ ^ g^^^ 
+4x(x + 2)n{n + l)Ho + 10x{x + 2){n + l)Ho 
+ {x-l){x-2)(2n+l)Ho . 

Therefore 

(2 + 3a; + 2nx)H2 - (2 + a; + 2nx)Hi2 - (8 + 16a; + 12na; + 3x'^ + 2nx^)Hi 

+ (2 + 4n + 17a; + 22nx + Sn^x + llx^ + l&nx^ + An^x'^)Ho = . (18) 

6.2 Equation from 7A([5'/(i + j + 1)]) 

Now apply 7 to the matrix obtained by evaluating the second identity (|15p at i + j + 1. If we denote 
this matrix by [SI{i + j + 1)], then this computation is the expansion of jAi[SI{i + J + 1)]) = from 

m- 

= a;7A([(i+j>j+j+3]) + (4a; + 2)7A([a,+j+3]) 

-x{x + 6)jA{[{i + j)a^+J+2]) - (4a;2 + 22x + 8)JA{[a^+J+2]) 
+4x{x + 2)jA{[{i + j)ai+j+i]) + Ux{x + 2)'yA{[ai+j+i]) 
+ {x - l)(a; - 2)7a([q+j+i]) - 4(x - l)(a; - 2)7a([c,+,]) . 

Using Tabled 

= x{2nH3-2{n-l)H2i+2{n^2)Hi3) + {4x + 2){H3- 

-x{x + 6){2nH2 ~ 2{n - l)i?i2) - {ix^ + 22x + 8){H2 - H12) 
+Ax{x + 2)2nHi + Ux{x + 2)Hi 

+ {x -l)ix- 2){2Hi + 2n(a; + 4)i7o) - 4(a; - l)(a; - 2)(2n + l)Ho ■ 

Therefore 

(1 + 2a; + nx)!!^ - (1 + a; + nx)H2i + (1 + nx)Hii 

-(4 + 11a; + Qnx + 2x^ + nx^)H2 + (4 + 5a; + Qnx + x^ + nx^)Hi2 (19) 
+ (2 + 11a; + 8nx + Sx^ + 4nx'^)Hi + (-2 + .x)(-l + x)(-2 + nx)Ho = . 
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6.3 Equation from 7A([ai+j+i], [SI{i + j)]) 

Now consider the expansion of 'yA{[ai+j+i], [SI{i + j)]) = from p^ . 

= x-fAi[ai+j+i], [(i +i)a,+j+2]) + (3.t + 2)JA{[a^+j+l], [0^+^+2]) 

-x{x + 6)7A([ai+i+i], [{i + j)ai+j+i]) - (Sx^ + 16x + 8)7A([a^+j+i], [a,;+j+i]) 

+4x(a; + 2)^Ai[ai+j+i], [(« + j)ai+j]) + 102:(a; + 2)-fAi[ai+j+i], [a^+j]) 

+ (x - l)(a: - 2)jA{[ai+j+i], [c,+j]) - 4(a: - l){x - 2)7A([a,+j+i], [q+j-i]) . 

Using the 7^([ai+j+i], *) computations from Table[31 we get 

= x{2nH2i~2{2n-3)Hi3) + {3x + 2){H2i-2Hi3) 
~x{x + 6)2(2n - l)Hi2 - (3a;^ + 16a; + 8)2Hi2 
+Ax{x + 2)n{n - \)Hi + l{)x{x + 2)nHi 

+ {x -l)ix- 2){{2n - l)Hi - {2n - l){x + 4)i?o) - 4(a; - l){x - 2)(-2ni7o)- 
Therefore for n > 2, 

{2 + 3x + 2nx)H2i - 4(1 + nx)Hi3 - 4(4 + 5.x + Qnx + x^ + nx^)Hi2 (20) 
+ (4n - 2 + 3x + Qnx + %n^x - x^ + ?>nx^ + An^x^)Hi - {x - 2){x - l)(2nx - 4 - x)Hq = 0. 

6.4 Two equations from the third identity 

The third identity is as given in Lemma [31 Define the column vector 

The third identity ([TB)) says that the vectors uq, f 1, . . . , Vn+2 are hnearly dependent with the weights 
given in (flT)) . i.e. 

n+2 

^ Wnj Wj = . (21) 

Now consider the determinant of the (n + 1) x (n + 1) matrix whose first n columns are the columns of 
A, and whose last column is the zero vector. Writing the zero vector in the form (PT|) and expanding 
the determinant by linearity, we find 

Substituting the weights from (|17p . this gives the equation 

(2 + 3x + 2nx)H2 - (10 + 4n + 15x + 16nx + An^x + 3x^ + 2nx^)Hi 

+ {n + l)(2n + 5)(2 + 2,x + 2nx + 2x^)Ho = . (22) 

Next we apply the same expansion trick to the matrix whose first n — \ columns are those of A, 
i.e. Wo, wi, . . . , whose {n — l)-st column is w„; and whose last column is the zero vector, written 

in the form (|2ip . Expanding the determinant by linearity, this time we obtain 

Wn,n+2H21 + W„,„+li/i2 - Wn,n-lHo = . 
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Therefore another equation is 

3(2 + 3.T + 2nx)H2i - 3 f 10 + 4n + 15a; + 16nx + 4n'^x + 3x'^ + 2nx^) H12 



(23) 



-(2n{l + 2n)(5 + 2n) + n{l + 2n)(3 + 2n)(5 + 2n)x + 3n(3 + 2n)(5 + 2n)x^^Ho = 



Equations (UHl), (HH), ([22]), form a 5 x 5 Hnear system Mu = b which expresses the 

determinants 

u = -ffs, H2i,Hi3 , H2, Hi2 
in terms of the two determinants Hq, Hi. The matrix M is as fohows: 





nx + 2x + 1 







with 





—nx — X — 1 
2nx + 3a; + 2 




2nx + 3x + 2 
2 





na; + 1 ~nx 

~4{nx + 1) 

2na; + 3a; + 2 



2x^ — 6nx — llx — 4 nx^ + x'^ + 6nx + 5a; + 4 



2na; — x — 2 

+ a^ + 6n; 
4 (na^ + a;^ + 6na; + 5a; + 4) 







3(2na; + 3a; + 2) 







-3(4a;n^ + 2a^n + 16an 
+ 4n + 3x^ + 15a + 10) 



det(A/) = 12(1 + na;)(l + 2x + nx){2 + x + 2nx){2 + 3.t + 2na;)^ . 
Solving Mu = b for u, we obtain each of iJa, H21, His, H2, H12 in terms of Hq and Hi. 



3(2 + 3a; + 2nx)H3 



-2(n+ 1) 8a;n 



12x^n^ 



64a;n 



8n^ + 6x^n 



+66x^71 + 162xn + 52n + 15x^ + 90x'^ + 126x + 84jiJo 
+3 (axti^ + Ax^n^ + 32^^^ + An^ + 2a;^n + ISx^n (24) 



+81a;n + 26n + 3x^ + IBa;^ + 63a; + 42 , 



3(2 + a; + 2nx)(2 + 3x + 2na;)i72i = ( 



64a;^ri'^ 



AlQx^n^ 



I2d>xn' 



192a;^n^ 



3(2 + a; + 2nx)H 
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(2 + 3x + 2nx)H2 
(2 + X + 2nx)Hi2 



-1040a;2n3 - 704x71^ - Un^ + I2x^n^ - 192x^n^ (25) 
-1192a;^n2 - 1360a;n2 - 288n^ + 24a;''n + 24x^n - 480x^n 
-1120a;n - 416n + 9x'^ + 63x^ + 48x'^ - 300a; - 240^ Hq 

+3(4a;n^ + 4xn + 4n - a; + 2) (^4a;7i^ + 2a;^n + 16a;n 

+4n + 3a;2 + 15a; + lO^iJi , 

(-16a;n'' - 32xn^ - 16n^ + 28x71^ - 2471^ - 6a;^77 - 12x^7i 
+80a;7i + 1671 - 3x^ - 12x^ + I2x + 48)Ho (26) 

+3 (4X77^ + 471^ + 2X^71 - 9X77 - 27t + x^ -4) Hi , 
-(77+l)(277 + 5) (2 x^ + 27ix + 3x + 2) Hq 

+(4xn^ + 2x^77 + 16x77 + 477 + 3x^ + 15x + 10)-ffi 



2x^77 — 9 X77 — IO77 + X^ 



(-4X77^ - 12X77^ - 477^ 
+2X ~ 8)Hq + (4X77^ + 4X77 + 477 - X + 2)Hl 



(27) 



(28) 



Equipped with these expansions, we now proceed with the calculation of the derivatives of Hq 
and Hi. 
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7 The derivatives of Hq and Hi 



7.1 The derivative of Hq 

From Definition [H 

Therefore by Proposition [2] we have 



Using FI{i + ]), 

{x - 2)x{x + 2)(3a; + 2)-iA{[FI{i + j)]) 



The values for 7/1 (*) from Table [2] give 

(x - 2)x{x + 2)(3x + 2)7^([F/(z + j)]) 



2(2: - 1)7a([(« + j)aj+i+2]) 

+ (a; - 6)(x - 2)7A([(i + j)aj+j+i]) 

+ (3a;2-2a; + 4)7A([ai+j+i]) 

-4(a;2 + 4)7^([(i+j>,+j]) 

-(3a;3 + 18a;2 - 2Qx + 24)7A([ai+j]) 

+8(a;-l)V([c.+j]) 
-32(a;-l)27^([c,+j_i]) . 



2{x - l){2nH2 - 2{n - l)Hi2) 

+ {x - 6){x - 2)2nHi 

+ {3x'^ ^2x + A)Hi 

-A{x^ + 4)n(n + l)Ho 

-{3x^ + 18x^ - 20x + 24) (71 + l)Ho 

+8{x-lf{2n+l)Ho . 



Now using the expressions in ((27|) and ([28| for H2 and _ffi2 in terms of Hi, Hq, we obtain 



where 



Q-T^o — QoHq + QiHi 
ax 



Q = {x~2){x + 2){2nx + x + 2){2nx + 3a; + 2) , 
Qo = -{n+l){lQx^n^ + -ix^n^ + A8x^v? +32x1? + 8x^n 

+mx^n + SOxn + 16n + 3a;^ + 12a;^ + 12x + ^ 

Qi = (2n + 3) {A,n^x^ + Anx^ + + 8na; + 4^ . 
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7.2 The derivative of Hi 

To differentiate Hi we use tfie expression 

Hi = 7A([a,+j+i]) 
from Table [21 From Proposition [2] we liave 



Tlierefore, to compute -^Hi 

JA{[a^+J+ll[Fn^+J)]) and 7^([F/(z + j + 1)]) 
are needed. For the first one of these 

(x - 2)xix + 2){3x + 2)7^([a,+,+i], [Flii+j)]) = 2{x ~ l)7^([a,+j+i], [(z + j)a,+,+2]) 

+ {x - 6){x - 2)jA{[ai+j+i], [{i + j)ai+j+i]) 
+ {3x^ -2x + A)jA{[ai+j+i], [ai+j+i]) 
~A{x^ + 4:)JAi[a^+J+l], [(z +i)a,+j]) 
-i3x^ + 18x^ - 20x + 2A)-fA{[a^+J+l], [a,+,]) 
+8{x - l)27A([a,+j+i], [ci+j]) 
-32{x - l f'yA{[ai+j+i], [ci+j-i]) . 

Using the entries in Table[3]for the "fA{[o-i+j+i], *) computations, we get 

{x - 2)x{x + 2)(3x + 2)7^([a,+j+i], [FI{t + j)]) = 2{x - l)(2niJ2i - 2(2n - 3)Hi3) 

+ {x - 6){x ~ 2)2{2n - l)Hi2 
+ {3x^ - 2a; + 4)2i7i2 
-4(a;2 + 4)n(n - l)Hi 
-{3x^ + 18x^ - 20a; + 24)nHi 
+8{x ~ l)^((2n - \)Bi - {2n - \){x + 4)i/o) 
-32(.T- l)2(-2niJo) . 

For the term 7^([F/(i + j + 1)]), we obtain 

(a;-2)a:(a; + 2)(3a: + 2)7A([F/(z + j + l)]) = 2(x - l)7A([(i + j)a,+j+3]) 

+2(2; - l)7A([ai+j+3]) 
+ (x - 6)(a; - 2)7A([(i + j)a.+j+2]) 
+2(8-5a; + 2a;2)7^([a,+j+2]) 
-4(x2+4)7A([(i + j)a,+j+i]) 
-(40 - 20a; + 22a;2 + 3a;3)7^([a,+j-+i]) 
+8(a;-l)V([c.+j+i]) 
-32(a;-l)27A([c,+,]) . 
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Using Table [2] this gives 

{x - 2)x{x + 2)(3x + 2)jA{[FI{t +j + 1)]) = 2{x - l)(2ni?3 - 2{n - l)H2i + 2{n - 2)H^,) 

+2{x ^ l)iH3 ~ H21 + H,s) 
+ {x ~6){x~ 2){2nH2 - 2{n ~ l)Hi2) 
+2{8-5x + 2x'^){H2- 
-4(a;2 + 4)2nHi 
-(40 - 20a; + 22x^ + 3x^)Hi 
+8{x - lf{2Hi + 2n{x + A)Ha) 
-?,2{x-lf{2n+l)Ho . 

Adding, we get 

[x - 2)x{x + 2){'ix + 2)-^Hi 
ax 

as a combination of H^, H21, Hi3, H2, H12, Hi, Hq. After that, we use the expressions (p4|) -(|28 |) for 
H3, H21, Hi3, H2 and H12 and express -^Hi as a hnear combination of iJp, Hi as 

U^Hi = UoHo + UiHi . (31) 



We find 



U = {x-2){x + 2){2nx + x + 2){2nx + 3a; + 2) , 
C/o = -2(n + 1) [iQx^n^ + Sx^n^ + 72x^n^ + 32a;n^ 

+2d,x^n^ + liea;^?!^ + 112a;n^ + IQn^ + 26a;^n + %Qx^n 
+104a;n + 56n + 7x^ + 22a;2 + 20a; + 56) , 



Ui = [iQx^n'^ + Ax^7i^ + Mx^TT' + 32a;n3 + I2x^n^ (32) 
+92x^n^ + 80a;n^ + I&t? + lla;^n + b&x^n + 44a;n + 32n 
+3a;^ + lOx^ - 4x + 24 

The exphcit polynomials in ([50]) and ([5^ are the coefficients of the system of differential equations 

(ED. 

Differentiating both sides of and substituting the expansions of ^i?o and ^-ffi in terms of 
Ha and Hi, we obtain 

R——^Ho = RqHo + RiHi , (33) 



where 



i? = (x - 2)2(a; + 2)^(2na; + a; + 2)(2nx + 3a; + 2) , 
Rq = {n + 1) (^4n-\T* + len^a;* + 19na;'' + 6.t'* + 32n^a;^ + 96n'^x^ 

+647ix^ + 18x^ - ASn^x'^ + 240nx^ + 48x'^ + 128n^x + 288n'^x 

+160na; + 264a; + 12871^ + 208n - 96^ , (34) 

Ri = -2{2n + 3) (^An^x^ + Anx^ + a;^ - An^x^ + 8na;^ - a;^ 

+16n^a; + 8nx + I2x + 16n - 4^ . 
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From dSni) and dMl), we find that QiR, RiQ, and RiQo ~ QiRo in 

QiR^Ho - RiQ^Ho + {RiQo - QiRo)Ho = (35) 

have GCD 

{2n + 3){x- 2){x + 2){2nx + x + 2){2nx + 3x + 2) . 

Dividmg through (j35p by this and defining S2, Si, Sq as the resuhing quotients, we obtain the second 
order differential equation satisfied by Hq. We record this in the following theorem. 

Theorem 2 Suppose the polynomials ak{x) and the (n + 1) x (n + 1) Hankel determinant Hq = 
H()[n,x) are as defined in ^ and Then 

S2^Ho + Si^Ho + SoHo = , (36) 
ax^ ax 

where 

S2 = {x-2){x + 2) {An^x^ + Anx^ + x^ + 8nx + 4) , 

Si ^ 2{Wx^ +Anx^ + x^ -Wx"^ + 8nx'^ -x"^ + l6n'^x + 8nx + 12x + 16n-A) , 
So = -n{n + 1) (4n^a;^ + 4nx^ + + 8nx - 8x + 36) . 

8 Evaluation at special points 

At this point we have enough information to evaluate Hq{x) at special points x without making use 
of the differential equation ([55)1 itself. 

Using the notation that incorporates the sizes of the matrices involved, we recall the following 
general result on Hankel determinants proved in [2]: 

Proposition 3 

Ho{n — 1, x)Ho{n + l,x) = Ho{n, x)H2{n, x) + Ho{n, x)Hi2 (n, x) — Hi{n, x)'^ . (37) 

8.1 Specialization at x = 2 

At X — 2, the derivative expression in gives 

~2{n + 1)(3 + 6n + 2n'^)Ho + (1 + 4n + 2n^)Hi = . 
From equations ((27)) and ([28l) at x = 2, 



{n + 2)H2 = -(ra+l)(n + 4)(2n + 5)i?o + (13 + lln + 2n2)iJi 
{n + l)Hi2 = -n{n + l)(2n + 5)Ho + n(2n + 3)Hi . 



Therefore at x = 2 we can write ((57|) as 

^2 - i)(7 + 8?i+ '^-2 
(1 + 4n + 2n2) 
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This is a recursion in i?o(n, 2)/Ho{n — 1, 2) with Hq{Q, 2) = 1, Ho{l, 2) = —7. Solving, we find 

Ho{n,2) = (-l)"(2n2 + 4n+l) . 
At a: = 2, the entries of the determinant in ([55]) speciaHze to 

a.(2) = 4'=+i-n' + 'V (38) 



The evaluation of the corresponding Hankel determinant is as follows: 
Corollary 1 Suppose ak{x) is as defined in Then 

Ho{n, 2) = det k+,(2)]o<,^^.<„ = (-l)"(2n2 + 4n + 1) . (39) 

8.2 Specialization at x = —2 

At a; = —2 the expression for the derivative in (l29l) gives 

-2(n + 1)(3 - 4n + 6n^ + An^)Ha + (2n + 3)(1 + 2n^)Hi = . 

Again from equations ([?7)l and we obtain at x — —2, 

(n+l)i/2 = -(n-l)(n+l)(2?i + 5)i7o + (n + 2)(2n+l)i?i , 
niJi2 = (2 - 4n - 571^ - 2n^)i7o + (n + l)(2n- l)i/i . 

Therefore we can use ([57]) at a: = — 2 and write 



(2n + l)(2n + 5)(3 - 4n + 2^2) (3 + 4n + 2n2) 
(2n + 3)2(1 + 2n2)2 



i/o(n - 1, -2)i/o(n + 1, -2) = ^ : - ^" ^ - > Ho{n, -2)2 



This is a recursion in Ho{n, —2) / HQ{n — 1,-2) with _ffo(0, — 2) = 1, i/o(l,— 2) = 5, which can be 
solved to give the simple product evaluation 

Ho{n,-2)^h2n + 3){l + 2n^) . (40) 
o 

Therefore 

Corollary 2 Suppose ak{x) is as defined in Then 

det [a,+, (-2)]o<,^^^.<„ - J(2n + 3)(1 + 2n2) . (41) 

The entries in (|4ip do not seem to have as simple an expression as the 0^(2) given in pS]) . although 
from the alternate expression for the generating function of the a^, we get the generating function 
of these numbers as 

1 2 



l-y-Ay^^yt 1 + ~ 2y(l + 4y) ' 

where t is the generating function of the Catalan numbers, as in the proof of Lemma [1] in Appendix 
11. 
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9 The differential equation solution 

Natural candidates for the expansion of the power series solution to the differential equation ([55]) 
are around x = 2 and x ~ —1. 

9.1 Solution at X = 2 

Putting 



Fo(a;)=^6fc(a;-2)'= , 
fe=o 

we find that the satisfy 

16fc(2fc+ l)(2n2 + 4n+ l)6fe 8(271^ + - lOfc^n^ + IS/cn^ - - IG/c^n + 26fcn 

-In - 3fc2 + 4fc - l)6fc_i + 2{%n'' + 20n^ - \%k^n^ 

+60fcn2 - 46n2 - IQk^n + 68fcn - 58n - Ak^ + 15fc - 14^fefe_2 
+ (n + 3 - fc)(fc + n - 2)(2n + \fh},--i 

for fc > 2 with 6^ = for fc < 0. From (|36p . we get 

n(n+l)(2n2 + 4n + 3) ^ 
''^ = 6(2n2+4n+l) ^° ' ^^^^ 

and therefore each 6^ is a multiple of 60 • It can then be proved by induction that 

_ 2^2 + 4n + 2fc2 + 1 /n + k\ 
^ " (2n2 + 4n+l)(2fc + l) 1^2^'° 

Since 60 = -ffo(2), 

i/o(2) ^ 2^2 + 4n + 2/s2 + 1 ^„ + fc^ 

fe=0 

The determinants at a; = 2 have the simple evaluation we already found in p9p . so that 

^.(.) . (-1). f; + + ' - 2)' ^ (43) 

The coefficients in (I43|) can be rewritten as binomial coefficients to obtain the expansion given in 
([7]) of Theorem [T] at x = 2. Note that the alternate notation Ho{x) in (subscript indicating the 
zero partition), is the (n + 1) x (n + 1) determinant denoted by Hn{x) in Theorem [TJ 

Using the expansion at x = 2, we can immediately write down the generating function of the 
Hn{x). We omit the proof of the following result. 

Theorem 3 Suppose ak{x) is as defined in Then 



^°(")=2n2+4n+l2:^ 2fcTl 2fc " 
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9.2 Solution at x = -2 

For the solution around x = — 2, put 

We find that the dk satisfy 



k=0 



I6k{2k + 3){2n^ + l)dk = - 8 (2n'^ + 2n^ - 10 k^n^ + 10 kri^ + 7 n'^ ~ ik'^n + 



6kn 

+5n- S/c-" + 2k + ljdk-i + 2[8n'' + 12n'' - lefc^n^ + 52kn^ 
-30^2 - Uk'^n + Ukn - Mn - Ak^ + 13fc - lo) rffc_2 
+ (fc - n - 3)(A: + - 2)(2n + lfdk-3. 



for fc > 2 with di. = for fc < 0. From 



we get 
n(l + n)(7 + 2n2) 



do 



10(1 + 2n2) 

and therefore each is a multiple of c?o. It can be proved by induction that 

+ k^ 



(45) 



du = {-lf 3(2^2 + 2fc2 + 4fc+l) 



(l + 2n2)(2fc + l)(2A; + 3) V 2fc 



do . 



Since do = Ho{-2), 



^ ^ Hai-2)^, ^,,3(2n2 + 2fc2 + 4fc + l) /n + fc 



l + 2n2^' ' (2/fc + l)(2fc + 3) 



2fc 



(a; + 2)* 



Using the evaluation of the determinants at a: = — 2 from (|40|) we obtain 

J, 2^2 + 2fc2 + 4/c + 1 + k 



iIo(a;) = (2n + 3)^(-l)* 



fc=0 



(2fc+ l)(2fc + 3) 



2fc 



{x + 2)^ . 



(46) 



The coefficients in p6)) can be rewritten in the form ([7]) of Theorem [T] Again, note that Hq{x) in 
is the (ri + 1) X (n + 1) determinant Hn{x) in Theorem[T] 
Evaluating ([7]) and ([8|) at x = we obtain the expressions 



det 



2(i+j)+2 
i+j 



E 

fc=0 



0<4j<n ;.^o 

n + fc 
2k 



[n + k + l) 



(2n + 3) 



(2n + 4fc + 1) 



n + k 
2k + 1 

n + k 
2k +1 



+ (2fc + 1) 



+ 8(fc + 1) 



n + k + l 
2k + 1 

n + k + l 
2k + 3 



{-2f (47) 
(-2)'= (48) 



which are alternate ways of writing the known evaluation of this determinant from ^ . 

As another corollary of Theorem [1] we have the following Hankel determinant evaluation at 
a; = 1, which depends on the residue class of n modulo 3: 



Corollary 3 



det 



2{i+j) + -i 
i+j 



0<i,j<n 



i(2n + 3) zfn = (mod 3) , 
-|(n + 2) ifn = 1 (mod 3) , 
i(2n + 5) ifn = 2 (mod 3) . 



(49) 
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Proof Since 



(50) 

the determinant is simply Ho{n, 1). We use the expression for the determinant ([7]) of Theorem [T] 
evaluated at a: = 1. Putting n = 3m, n = Sm + 1 and n = 3m + 2 for the three residue classes 
modulo 3, the Corollary is a consequence of the resulting binomial identities 



3m 



2m + 1 = (-1)"^ 
fe=o '- 

3m+l 

4(m + l) = i-ir E 

fe=0 

3m+: 

2m + 3 = (-1)" J2 



^^"^ + ^)l2fc + l +(2^ + 1) 2fc+l 



fe=0 



, , 3m + fc + 1 , 
(6™ + 5)( + i )+(2fc + l) 

, , 3m + fc + 2 , 
(6™ + 7)( + i )+(2fc + l) 



(-1)' 

3m + fc + 2 
2fc + 1 

3m + fc + 3 



(51) 



(-1)' (52) 



2fc + 1 



(-1)' (53) 



which can be proved by making use of the generating function given in Theorem [3] at x = \. 



9.3 Solution at x = 

The power series solution to p6p around x is more difficult to derive directly. For and the 
Hn{x) as defined in ([5]) and ([6]) this expansion is given by 



H„{x) = ^(_l)"(»-l)/2+Mfc-l)/2+fcn (2fc + (-1)"-^-) 



L^JIfc! 



(54) 



We are grateful to the anonymous referee for pointing out the above explicit form of the determinant 
around x — 0. This expansion is an immediate consequence of the generating function for the 
determinants at arbitrary x that we have provided in (I44|) . 

Following the route of the proofs of the cases x = 2 and x = —2, one would put 



Hn{x) = 

and show that the satisfy the recursion 



ekx'' , 



k=0 



16fc(fc-l)efe = -8(fc- l)(4fcn- 12n+ l)efc-i 

-4 (4n^fc^ + 4nfc^ - 28n^fc - 24nfc - 6fc + 49n^ + 41n + 12) ek-2 
~2 (4n3 + 4fcn2 - Un"^ - Ak'^n + 20fcn - 28n + fc - 3) efc^g 
+ (fc - n - 4)(fc + n - 3)(2n + l)^efe_4 

for k > 2 with = for fc < 0. In this case each Cfc is a function of eo and ei. We know eo explicitly 
by However in this case a relationship similar to (|^^ and (|^5)) of the a; = 2 and .t = — 2 cases 
does not drop out of the differential equation to give a similar relation between bq and ei. This is 
because the special values of x that kills off the second derivative term in Theorem [2] are x — ±2. 
An alternate approach is to show directly that the coefficient of x'^ in (IM)) satisfies the recurrence 
for the Cfc , but again this would fall back on the already proved expansions of Theorem [T] for the 
value of the derivative at a; = 0. 
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10 Zeros of HQ{n,x) 



The determinants _ffo(n, a;) of Theorem [T] are not orthogonal polynomials. But they satisfy a re- 
currence relation with polynomial coefficients involving three consecutive terms of the sequence as 
follows: 

Corollary 4 

(2 + (2n + 3)xfHQ{n + 2,x)+ x(4 + 4(2n + 3)x + (2n + 3)(2n + 5)x^)HQ{n + l,x) 

+ {2 + {2ji + 5)x)^Hq{ji,x) =0 . (55) 

Proof The recurrence relation can be verified by making use of the explicit form of Ho{n, x) from 
Theorem [TJ • 

Table [T] gives a list of the zeros of Ho{l,x) through Ho{7,x). The zeros are real and interlacing. 
It is possible that the polynomials HQ{n, x) can be obtained from an orthogonal family by a suitable 
transformation. 

-0.333 
-0.358 0.558 
-0.601 -0.194 1.224 
-1.083 -0.207 0.324 1.522 
-1.367 -0.351 -0.137 0.815 1.678 
-1.540 -0.746 -0.146 0.229 1.127 1.768 
-1.651 -1.028 -0.246 -0.107 0.608 1.333 1.825 



Table 1: Zeros of the Hankel determinants Ho{l,x) through Ho{7,x) of Theorem[TJ 

A sequence of polynomials {Pn{x)}n>Q with degP„ = n is called a Sturm sequence on an open 
interval (a, b) if P„ has exactly n simple real zeros in (a, b), and for every n > 1, zeros of Pn{x) and 
Pn+i{x) strictly interlace. 

Theorem 4 Suppose and the Ho{n,x) are as defined in ^ and (0). Then {Ho{n, x)}n>o is a 
Sturm sequence on (—2,2). 

Proof 

Consider the two expansions of HQ{n,x) in (|43|) and (|46)l . The first one of these implies that 
(— l)"iJo(n, x) > for X > 2, and the second one implies that HQ{n,x) > for a; < —2. Therefore 
the zeros of Ho{n, x) are contained in (—2, 2). 

We next prove that like orthogonal polynomials, HQ{n, x) has n distinct real zeros and the zeros 
of H(){n, x) lie strictly between the zeros of Hq{ji + l,x). This interlacing property is a consequence 
of the form of the recursion (|55p 

a'^Ho{n + 2,x)+xl3Ho{n + l,x)+j^Ho{n,x) = (56) 
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where /3 > for every x and n. We use induction on n. For any two consecutive zeros ri,r2 
of ffo(" + l,a;) the induction hypothesis impHes that Ho{n,ri) and Ho{n,r2) have opposite signs. 
Therefore from the recursion, HQ{n + 2, ri) and Ho{n + 2, r2) also have opposite signs and so Ho{n + 
2,x) has at least one zero in the interval (ri,r2). This accounts for > n zeros of Ho{n + 2,x). Let 
(52 < 2 be the largest zero of Ho{n + l,x). By the induction hypothesis, Ho{n,x) has no zeros on 
[62, 00). Therefore its sign at x ~ S2 is the same as its sign at a: = 2, which is (—1)". But the sign 
of Ho{n + 2, x) is also (—1)" at x = 2, but opposite of the sign of Ho{n, x) at x = S2 by ([SS]) . This 
forces Hfj^n + 2, x) to change sign and have a zero in (^2, 2). By a counting argument, Ha{n + 2, x) 
has to have another zero in (—2, 5i) where 5i is the smallest zero of i?o(f^ + 1, x). • 

11 Discussion, patterns and conjectures 

We introduced a class of multilinear operators 7 acting on tuples of matrices to take the place of the 
trace method of our earlier calculations. This approach to evaluate Hankel determinants is easier 
to work with: the 7-operators are easier to differentiate, and they do not produce the extraneous 
nonlinear terms. In the (2, 2)-case that we have covered in detail, we have also obtained numerical 
evaluations at special points as a byproduct. Furthermore we saw that the resulting polynomials 
have intriguing properties. 

Even though the application of the 7-operator reduces the calculations involved in almost product 
evaluations of Hankel determinants considerably, there are still stumbling blocks in the general 
(2,r)-case, and other cases that differ little from this. We consider a few of these determinants and 
conjecture closed forms for the evaluations. 

Corollary [3] is just one example of a strange pattern that holds for Hankel determinants where 
the entries are the polynomials a"^'^\x) defined in Taking a; = 0, let 




parametrized by r > 0. For notational simplicity, define 



F{n,r) = det [at+j] 



-no<i.,j<n ■ 



Then the evaluation in Corollary [3] can be written as 



F(3m,3) 
F(3to+ 1,3) 
F(3to + 2,3) 



2m + 1 



2m + 3 . 



4(m + 1) 



As an example, consider the following evaluations for the case r 



7: 



F(7to + 1,7) 
F(7to + 2,7) 
F(7m + 3, 7) 



F(7m,7) 



(2m + 1)3 

{m + l)(2m + l)2(9604m3 + 9604m2 - 1323m, - 2340)/90 
-(m + 1)2 (2m + l)(19208m3 + 67228m2 + 70854m + 23445)/45 
64(m + 1)3 
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F(7m + 4,7) = (m + l)2(2m + 3)(19208m3 + 48020TO2 + 32438m + 3015m)/45 
F(7m + 5,7) = -(m + l)(2m + 3)2(9604m^ + 48020m2 + 75509m + 38110)/90 
F(7to + 6,7) = (2m + 3)3. 

These evaluations have been verified for a significant range of m. This unusual set of formulas 
is typical of a complex pattern of evaluations of F{n,r) that continues with several unexpected 
dependencies on the value of n modulo r and on r modulo 4. For example, if r is odd then there is 
strong experimental evidence that 

F{rm, r) = F{rm - 1, r) = (2m + _ 

When we consider even r there is another twist to take into account. Experimental evidence tells 
us that 

[ 1 if r = (mod 4) , 

F{rm, r) — F{rm — 1, r) = < 

[ (-1)" if r EE 2 (mod 4) . 

Another interesting pattern we observe is the following for odd r: 

F{rm + (r - l)/2, r) = 2'^-i(m + l)('^-i)/2 . 

For even r there is also a simple pattern of this type: 

(-l)''/4+i(2r(m + ifr = (mod 4) , 



F{rm + r/2, r) 



(_l)('-+2)/4+m(-2^(^4_ l))r/2-i ^ ^ 2 (mod 4) 



In addition to these nice evaluations there are many that are not so simple. For example the 
F{rm + 1, r) becomes more and more complex as r increases. For r = 5 

F(5m+ 1,5) = -(m+ l)(2m+ l)(50m + 39)/3. 

For r — 7 the evaluation contains a cubic factor: 

F(7m +1,7) = (m + l)(2m + lf{960im^ + 9604m2 - 1323m - 2340)/90 

and when r — 9 the evaluation contains a quartic factor: 

i^(9m, + l,9) = -(m+ l)(2m+ I)3(3m + 2)(52488m^ + 69984m3 + 22518m2 + 1674m + 1505)/70 . 

We suspect that this irreducible factor keeps gaining a degree when r is increased by 2. 
These conjectures appear to be difficult to prove in their full generality using either the methods 
described in Krattenthaler [6l [7| or with the methods of the present paper. For any fixed r, the 
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methods of this paper might apply but it is hard to see how to approach the problem when r is left 
as a parameter. 

Further experimental evidence suggests that the determinants 



satisfy second order differential equations. However as r gets larger the differential equations and 
the first and second identities of our method become increasingly complex. We mention that there 
are also difficulties in evaluating the family of determinants 



For this family, the order of the differential equation for the determinant seems to increase with r. 
When r = 4, for example, experiments suggest that ([57)) satisfies a fourth order differential equation. 
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12 Appendix I 

The results in this Appendix apply to general Hankel matrices. We let x(S') denote the indicator of 
the statement S: x{S) = 1 if 5 is true and = if 5 is false. 

12.1 Properties of the 7-operator 

Proposition [1] 

For m < n + I, 

7^(Xi,...,X„) = ^det(As,.) (58) 

where the summation is over all subsets S of {0, 1, . . . ,rt} with \S\ = m and all permutation a of 
{l,2,...,m}. 

Proof Expand 

det(^ + tiXi + <2^2 + • • • + im^m) 

by columns (or rows) using the linearity of the determinant to obtain 

det(^ + tiXi + t2X2 + ■■■ + tmXm) = ^1*2 ' ' ' X! det(As,<,) (59) 

where As^cr is as defined in Definition [2l The proof follows by applying dt^dt2 ■ ■ ■ dt^ and putting 
ii = • • • = t„, =0. • 

Proposition [2] 

For m < n, 

d . \ / di . , ^ , d . 

'J^^a{Xi, . . . ,Xm) = jA[-^A,Xl, . . . ,Xm) + 2_^"fA{Xi, . . . ,Xj-i, —Xj,Xj + i, . . . , X^) 

Proof By Proposition [T] and the expression in , 

S.(T 



Let B = As^a- By Cramer's rule, 

d 1 " 

T^(^"^^) = d^g^^*(^^) 

where Bj is obtained from B by replacing the j-th column of B by its derivative. In terms of the 
matrix A, let As^aj denote this matrix. 
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Therefore 



^^det(As,.,,) 

n n 

E E ^ det(As,.,,) + E E e ^) det(As,.,,) 

j=0 S,cr j=0 S,a 

d d 



12.2 Expansion of the convolution matrices 

The expansion of the convolution matrices [ci+j+fe] for /c > — 1 are as foUows: 
Proposition 4 Suppose the convolution polynomial Cn is as defined in Definition\^ Then 

n+k 

[Ci+j+k\Q<i,j<n = Qp[Qi+j+fc-pX(j ^ P ~ fc)]o<»,j<K 



(60) 



p=0 



n-l 



p=0 

Proof The (i, j)-th entry of the matrix on the right-hand side of (I5D1) is 

n+k n— 1 

The upper hmit of the sums need not go past i + j + fc. In the second sum, replace phy i+ j + k — p 
and rearrange the indices. We get 

i+j+k i+j+k i+j+k 

E o,pai+j+k-pX{j >p-k)+ E o-pO-i+j+k-pXij <p - k) = E apO-i+^+k-p 

p—0 p—0 p—0 



Below are a few examples of the expansion of the convolution matrices obtained from (|60p . For 
k = -1, 








ao 


ai 










ao 


Ci+j-l]o<iJ<2 — flO 





ai 




+ fli 








ai 







0.2 


"3 . 










02 





ao ai 02 
ai 02 03 



Ol 






Oq Ol 02 



28 



For fc = 0, 



'i+j]o<iJ<2 — OO 



ao ai 02 
ai 02 03 + 0-1 
02 03 04 



flO Oi 02 03 

O2 O3 04 



oo oi 
oi 02 
02 03 



a2 



00 
01 
02 



Oi 







01 02 03 



For A: = 1, the expansion is as given in ([Tc 
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13 Appendix II 

In this Appendix we give the proofs of the statements needed for the (2, 2)-case. 

The proofs of Lemma [TJ Lemma [5] are based on generating function manipulations, as given 
below. The first identity for the (2, 2)-case is: 
Lemma [T] 

(x - 2)x(x + 2){3x + 2)4^a„ + (3x^ + ISx^ - 20a; + 24 + 4n (x'^ + 4)) a„ 
dx 

- {n{x - 6)(x - 2) + 3x'^ - 2x + A) a„+i - 2n{x - l)a„+2 (61) 

-8{X - ifCn + 32{X - lfCn-1 = 

Proof From [121 E] , we have 

+3 

/ = f{x,y) ^ak{x)y 

Here 

fc>0 ^ ' 
satisfies 

yt^ = t-l . (62) 
Using -j-t — I {\ — 2yt) in the computation oi and using the resulting expressions for and 



{2-m-xyt) 



/' = 4-f, we make the substitutions 



d 

dx 


dx 




an - 


- f 




nan - 


- yf 






- (/- 




na„+i - 






nan+2 - 




-l-(4 


^n ~ 













,2\' 



in the left-hand side of (pT|) . The resulting expression factors as 

{t-i-yt') 



(t-2)V(l-2ty)(l-tey) 
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(642:2?;3t^ - 128xyh^ + 64y^t^ - IQx'^y'^t' + 32xyh^ - lOy^i^ 



+20x^y^t^ - IGx'^yh'^ + 16xyh^ ~ x^yH'^ - ix'^yh^ ~ 12xyh^ + 32yh^ + Gx'^yt^ 
+2xyt'^ ~ 8yf^ - 32x^y^t^ + 32x^y^t^ - Gx^y'^t^ - ix'^y'^t^ + iOxy'^t^ - 80y^t^ 
-Axt^ - Ax'^yt^ - Axyt^ + 8yt^ + At^ + l&x'^y'^t^ + A&x^y'^t^ - GAxy^f + 4xf 
-IQxyt^ + l&yt^ - At^ - 32x'^yt + 32yt + I6x - 16^ 

and therefore vanishes by ([62|) . • 

The second identity is: 
Lemma [2] 



{nx + 3x + 2)a„+2 — {tix{x + 6) + 3x^ + 16x + 8)a,i+i + 2x{x + 2){2n + 5)a„ 

+{x -l){x- 2)c„ - 4:{x - l)(x - 2)c„_i = (63) 

Proof Again passing to the generating functions, we find that the generating function of the 
left-hand side of (l63l) factors as 



{t 1 yt^) foA^n,3-t5 o^2„,3.5 iR„,3.5 , o^2„,2.5 



(24:xyh^ - Sx^y^t^ - 16y^i^ + 2x^y^t^ - 6xyh^ + Ayh^ 



(i_2)V(l-2ty)(l-tey)2 
+8x^y^t'^ + Wx'^y^t^ ~ x^y'^t^ ~ Gx'^y'^t^ ~ Axy'^t'^ ~ 8y'^f^ + xyt^ + 2yt'^ - 8x^y^t^ 

-Wx'^y^t^ - Ux^y^t^ - Wxy'^t^ + 16y^t^ - xf + 2x^yt^ + lOxyt^ + Ayt^ - 2t^ + Ax^y'^t^ 

+24:x'^y^f + 32xy^t^ + xt^ - Axyt^ - 8yt^ + 2t^ - 8x^yt - 2Axyt - IQyt + Ax - 

which again vanishes by ([5^ . 

The third identity is: 
Lemma [3] 



n+2 

w« J {x)ai+j {x) = (64) 

3=0 



2{2n + 5) /n + j + 2\ (2n + 3)(2?i + 5) /n + j + 2 



for i = 0, 1, . . . , n where 

WnA^) = (-1)' 



2j + 3 V 2j + l 

We do not give the proof of the third identity Lemma [3] but remark that once the weights are 
guessed, the proofs of the identities can be left to automatic binomial identity provers such as 
MultiZeilberger supplied by Doron Zeilberger (in Maple jl4j), and MultiSum by Wegschaider (in 



2j + l V 2j ) 2j + l V 2j 

^(2n + 3)(2n + 5)/n + , + 2X ,| ^^^^ 
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Mathematica [13]). The main step in finding the coefficients is interpolation and a symbolic algebra 
system (Mathematica in our case). 

The weights in general can be found from the relation 

Ho = (66) 

which holds for k — 0, 1, . . . ,n. This can be seen by computing the determinant of the matrix 
obtained from A — [cLi+j]Q<i.j<n by replacing column n — fc by column n, and column n by the zero 
vector written as sum of column vectors as indicated by the third identity. In the present case this 
is (j64p . Expanding, all but three determinants vanish, giving (|66p . 

We use (j66| to guess third identities in general. For instance with offset 2 (i.e. the vectors 
involved in the third identity are vq through Vn+2), it is possible to first guess Wn,n+2,Wn,n+i,Wn,n 
by linear algebra, then use (j66p to solve for Wn^n-k and consequently find the candidate coefficients 
by interpolation. 
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14 Appendix III: Tables of 7-operator evaluations 



The tables given in this Appendix apply to general Hankel matrices. 

lAilai+j]) = {n + l)Ho 

7A([ai+i+i]) = Hi 

7A([ai+j+2]) = H2 — H12 

lA{[ai+j+z]) = H3-H21+H13 

7A{[ai+j+4]) = H4- Hsi + H212 - Hn 

lA{[ai+j+b\) = Hr-, - H41 + H^i2 — H213 + Hia 

lA{[{i + i)ai+j]) = n{n+l)Ho 

lA{[{i + 3)ai+j+i\) = 2nHi 

-lA{[{i+3)ai+j+2]) = 2nH2 - 2{n - l)Hi2 

7A([(i + 3)ai+j+a\) = 2nH3 - 2(n - l)H2i + 2(n - 2)Hi3 

7^([(i + j)a,+j-+4]) = 2nHi - 2{n - l)H3i + 2{n - 2)H2i2 - 2(n - 3)i?i4 

7A([(i + J>^+J+5]) = 2n7?5 - 2(n - l)i?4i + 2{n - 2)H3i2 - 2(n - 3)iJ2i3 + 2(n - 4)Hi5 

lA{[Ci+j-l]) = 

lAiici+j]) = {2n + l)aoHo 

7A([ci+j+i]) = 200^1+27101^0 

jA{[ci+j+2]) = 2aoH2 - 2aQHi2 + 2aiHi + {2n - l)a2Ho 

7A([cj+j+3]) = 2aoH-3 - 2anH2i + 2aoHi;i + 2aiH2 
-2aiHi2 + 2a2Hi + (2n - 2)03^0 

Table 2: 7a(*) computations. 
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lA{[ai+j+i], [fli+j+i]) 
7A([ai+j+i], [ai+j+2]) 

lA{[a^+j+l],[{i + j)a^+j]) 
7A{[ai+j+i], [{i + j)ai+j+i]) 
lAiiai+j+i], [{i + j)a,+j+2]) 
7A{[ai+j+i\, [{i + j)ai+j+3]) 

lAi[a^+j+l], [{i + j)a^+j+i]) 

JA{[a^+j+l], [Ci+j-l]) 

lA{[ai+j+i], [c,+j]) 
lA{[ai+j+i], [ct+j+i]) 

lA{[ai+j+l],[c^+j+2]) 



= nHi 

= 2Hi2 

= H21 — 2Hi?, 

— H31 — H22 — + 1H\i 

— Hm — H32 — H312 + H221 + H213 — 2Hi5 
= n{n - l)Hi 

= 2(2n - l)Hi2 

= 2nH2i - 2{2n - 3)Hi3 

= 2ni?3i - 2(n - l)i?22 

-2(n - l)i?2i2 + 2(2n - 5)i?i4 

= 2nH4i - 2{n - l)i?32 - 2(n - l)i?3i2 

+2(n - 2)if2n + 2(n - 2)^/213 - 2(2n - 7)Hi5 

= — 2nao-ffo 

= (2n - l)aoi/i - (2n - l)aii?o 

= 4aoifi2 + 2{n - l)aiHi - 2{n - l)a2Ho 

= 2aoi?2i - 4aoi?i3 + 4aiHi2 + (2n - 3)02-^1 - (2n - 3)a3i?o 

Table 3: ^A{[ai+j+i],*) computations. 



7A([ai+j+2],[ai+j-]) 
7A([ai+j+2], [aj+j+i]) 
7a( [ai+j+2], [aj+i+2]) 
7A([ai+j+2], [aj+j+a]) 
7A([ai+j+2], [(z + j)a^+j]) 
7A([ai+j+2], [(« + j)ai+j+i\) 
7A{[ai+j+2], [{i + j)ai+j+2]) 
7A{[ai+j+2], [{i + j)ai+j+3\) 
7A([a^+j+2], [ci+j-i]) 

lA{[ai+J+2],[c^+j]) 



= nH2 — nH 12 

= H21 — 2Hi3 

= 2H22 — 2H212 + 2Hii 

= H32 - H312 -H221+ 2H213 - 2Fl5 

= n(n - 1)H2 - (n^ - n + 2)Hi2 

= 2{n - l)H2i - 4(n - l)Hi3 

= 2{2n - 1)H22 - 2{2n - 2)^^212 + 2(2n - 4)ifi4 

= 2ni?32 - 2ni?3i2 - 2(n - 2)iJ2n + 4(n - 2)iJ2i3 - 4(n - 3)i?i5 

= -2aoi?i - 2{n ~ l)aii?o 

(2n - l)aoH2 - (2n - l)aoHi2 - 2aiHi - (2n - 3)a2i/o 

Table 4: 7A([aj+j+2], *) computations. 



^A{[ai+j+i], [aj+j+i], [ai+j]) 
7A([ai+j+i], [a^+j+i], [aj+j+i]) 
7A([ai+j+i], [oii+j+i], [ai+j+2]) 
7A([ai+j+i], [a^+j+i], [ai+j+3]) 
jA{[ai+j+i], [ai+j+i], [(z + j)a,+j]) 
lAiiai+j+i], [ai+j+i], [{i + j)ai+j+i]) 
7Ai[ai+j+i], [aj+j+i], [{i + j)aj+j+2]) 
nfA{[ai+j+i], [ai+j+i], [(i + j)ai+j+3]) 
jA{[at+j+i], [ai+j+i], [ci+j-i]) 
7A([ai+j+i], [tti+j+i], [ci+j]) 



2{n - l)Hi2 
6Hi3 

2H212 — QHii 

2H312 — 2_ff22i — 2iJ2i3 + 6Hi5 

2(n-l)(n-2)7?i2 

12(n - l)i?i3 

AnH2i2 - 12(n - 2)Hi4 

AnH3i2 - 4(n - l)iJ22i - 4(n - I)il2i3 + 12(n - 3)i?i5 

-4(n - l)ao-fri + 4(n - l)aii?o 

2ao(2n - 3)i/i2 - 2(2n - 3)aiHi + 2{2n - 3)a2-ffo 



Table 5: 7A([ai+j+i], [oj+j+i],*) computations. 
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